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The modification of acoustic phonons in semiconductor nanostructures embedded in a host crystal 
is investigated including corrections due to strain within continuum elasticity theory. Effective elastic 
constants are calculated employing ab initio density functional theory. For a spherical InAs quantum 
dot embedded in GaAs barrier material, the electron-phonon coupling is calculated. Its strength is 
shown to be suppressed compared to the assumption of bulk phonons. 
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I. INTRODUCTION 

Semiconductor nanostructures have attracted increas- 
ing interest over the last couple of years. Especially quan- 
tum dots (QD) are proposed as possible candidates in 
applications like single photon sources and quantum in- 
formation devices. Improved growth techniques nowa- 
days enable the fabrication of small QDs with a narrow 
size distribution. Single, isolated QDs embedded in a 
semiconductor matrix or in solution are of high inter- 
est for basic research. Since their electronic properties 
are inevitable connected to the underlying lattice and its 
vibrations (phonons) they offer the possibility to study 
these fundamental interactions. While the theory for 
electron-phonon interaction in bulk systems is well es- 
tablished, the inhomogeneous nature of nanostructures 
leads to strong modifications of the electronic properties 
as well as the phonon spectrum. Usually, QDs are de- 
signed such that at least the energetically lowest state 
is spatially confined. The existence of boundaries be- 
tween the constituting materials and/or to the vacuum 
introduces a coupling of the longitudinal and transverse 
phonon modes even for isotropic media. Additionally, 
new types of confined interface and surfaces modes can 
occur. Especially the electron-acoustic phonon interac- 
tion is of high interest since it is the source of so-called 
pure dephasing of optical excitations in QDs, as detailed 
in recent articles^. 

The electron-acoustic phonon interaction is usually 
treated by deformation potential coupling, and acoustic 
phonons are described as bulk phonons. Only recently 
some progress has been made to include phonons in in- 
homogeneous media in the treatment of electron-phonon 
interaction in slab and half space geometries 1 . So far to 
the best of our knowledge no investigations exist treating 
the modified phonons of the nanostructure itself in the 
interaction. 

The present article provides insights into the role of 
the phonon modification due to the nanostructure, in- 
cluding strain. We restrict ourselves to the simplest case 
- a spherical QD - in order to describe the qualitative 
effect on the electron-acoustic phonon coupling. The de- 
tails of the model and the corresponding theory are devel- 
oped in Sec.[n] Calculations of strain dependent acoustic 



phonons on the basis of ab initio density functional the- 
ory (DFT) calculations are given in Sec. lIIBl which serve 
as input for the calculation of the electron-phonon cou- 
pling function. Results are given in Sec. IIIII 

II. MODEL AND THEORY 

Epitaxial growth methods enable the fabrication of 
semiconductor heterostructures with sizes of a few 
nanometers only. The different materials are connected 
pseudomorphically, i.e. without structural defects. De- 
fects would lead to partial relaxation of intrinsic strain 
which results from different lattice constants of the con- 
stituents. The strain itself causes the change of physi- 
cal properties compared to relaxed bulk materials. This 
includes the electronic (band structure) as well as the 
lattice properties (phonons). 

The system we study in this article is a spherical shape 
inclusion (quantum dot) in an infinite matrix (barrier). 
This artificial system should serve as test model for the 
relevance of non-bulk phonon modes in the electron- 
phonon coupling. For the material parameters we choose 
the prototype InAs/GaAs system. The first part is de- 
voted to the description of the acoustic phonons within 
continuum elasticity theory followed by the description 
of the electron-acoustic phonon coupling. 

A. Acoustic phonons within continuum elasticity 

In elastic continuum theory, acoustic phonons follow 
from the wave equation of the displacement field u(r) 
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here written in Cartesian coordinates j, k — x,y, z. Both 
the mass density p(r) and the stress tensor <j(r) are in 
general spatially dependent. In linear elasticity theory 
the stress tensor a is proportional to the strain tensor 
e defined by the spatial derivatives of the displacement 
fields, both quantities being related by the fourth order 

tensor of elastic constants C (Hooke's law). For isotropic 
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systems there are only two independent elastic constants 
Cn and C44, which are related to the well-known Lame 
coefficients C12 = C\\ — 2C44 = A and C44 = /1. The 
stress tensor elements are then given by 
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The restriction to isotropic media allows to treat the 
phonons in large part analytically. 

For our spherical model system it is convenient to 
transform to spherical coordinates. For the specific 
electron-phonon coupling considered below, the only rel- 
evant phonon mode has angular momentum zero (breath- 
ing mode), u(r) = e r u(r), which fulfills the equation 
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with elastic constants and mass density having an arbi- 
trary radial dependence. This one-dimensional second- 
order differential equation contains the boundary condi- 
tions for a stepwise constant medium as well, which are 



u and C\i 



du 
dr 



2u 
r 



4C44 — continuous .(4) 



We solve Eq. @ numerically, starting with a regular so- 
lution u(r) oc r at the origin. The general normalization 
condition of the displacement within a large volume f2 
reads 3 



d 3 rp(r)u„(r) • iv(r) = S v y- — , 



(5) 



where v is the discrete mode index. This normalization 
is converted into a matching condition at large distances 
with an outgoing wave containing a phase shift, as stan- 
dard in scattering theory. As a consequence, v is replaced 
by a (continuous) energy variable tuo v — » E. 

The elastic constants in Eq. © are not to be taken 
as the usual tabulated bulk values for the dot and the 
barrier material: The phonon calculation has to start 
from the reference state where the crystal is locally under 
finite strain. Therefore the spatial dependency is not only 
due to the change in material but also due to the spatially 
varying strain. 



B. Effective Elastic Constants from DFT 

The usual treatment of continuum elasticity assumes 
fixed elastic constants. In our case, however, the system 
contains spatial regions where the material is strained. 
Intuitively, one would expect that a material under com- 
pression would become harder due to the presence of an- 
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FIG. 1: Effective elastic constants for GaAs (black) and InAs 
(red) in dependence on the lattice constant a for hydrostatic 
deformations calculated by DFT employing perturbation the- 
ory. The (experimental) equilibrium lattice constants are 
marked by full vertical lines, while the lattice constant of the 
relaxed InAs QD is given by the dashed line. 



harmonic terms in the expansion of the free energy den- 
sity around the equilibrium position, 

1 6 1 6 
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The equilibrium density is written as po- We restrict 
ourselves to hydrostatic compressions in the following. 
An extension to general deformations is straightforward 
but numerically a far more complex task. Allowing only 
small distorsions of the lattice (an expansion up to second 
order is sufficient), only the bulk modulus B — l/3(Cn + 
2C12) would be a constant quantity. Including up to 
third order in the expansion an analytical derivation gives 
a linear dependence on the strength of the hydrostatic 
lattice distortion^. 

The calculation of elastic constants is possible directly 
by employing DFT. Starting from a relaxed unit cell 
under the application of finite deformations within den- 
sity functional perturbation theory one can calculate the 
fourth order tensor of elastic constants. The density func- 
tional calculations are carried out with the ABINIT com- 
puter cod o 5 i 6 i 7 . The local density approximation is ap- 
plied for the exchange-correlation. Soft norm-conserving 
pseudopotentials are taken from the code of the Fritz- 
Haber Institute, Berlin^. Wave functions are expanded 
into plane waves with converged Monkhorst-Pack meshes 
offc = 8x8x8 per lxl unit cell and a cutoff energy of 
E cut = 36 Ha = 979 eV. 

The dependence of the calculated elastic constants on 
the (variable) lattice constant, which corresponds to hy- 
drostatic deformation, is displayed in Fig. ^ ln the range 
from -7% to +7% lattice constant variation we find a 
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nearly linear reduction of Cxi and C12 with increasing 
the lattice constant, which is expected from a third or- 
der expansion^, showing that even for small distorsions 
the assumption of constant elasticity does not hold. The 
deviation from the linear behavior, especially for C44, is 
a clear indication that even higher orders than the third 
one influence significantly the effective elastic constants. 
Numerical deviations from the experimental values are 
mainly due to the known underestimation of the lattice 
constant in LDA for III-V semiconductors (see Tab. QJ. 
Given a local static lattice distortion, it is now possible to 
implement these elastic constants in the phonon calcula- 
tion, which depend on the spatial position due to change 
in material and strain as well. 
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C. Electron-acoustic phonon coupling 

The electron-phonon interaction has two main con- 
tributions: Deformation potential coupling and piezo- 
electric coupling. The second contribution will not be 
considered here, although some of the following discus- 
sions apply also to this part. 

For conduction electrons at the T-point in III-V semi- 
conductors the deformation potential couples only to a 
volume deformation (divu). The interaction Hamilto- 



u„(r) + h.c. 



(6) 



with the (material-dependent) deformation potential 
constant for the conduction band Dc and the phonon 
creation operator b\,. In a QD with confinement func- 
tions ip n (r) of the electron, the relevant electron-phonon 
matrix elements are then 



M': 



d 3 r<(r)^ m (r)^ c (r)divu„(r). (7) 



The spherical symmetry of the QD leads to well-defined 
angular quantum numbers of the confinement levels, 
which select an appropriate symmetry of the lattice dis- 
placement. For electronic s-s-transitions, only the az- 
imuthal quantum number I — is needed (breathing 
mode). For the polaronic modification of a single elec- 
tronic transition within the independent Boson model 9 , 
the diagonal element in Eq. (0 is of central importance. 
Concentrating on the lowest (s-type) confinement state 
n = 1, we need to calculate the following coupling func- 
tion^ 



f{E) = Y,5{E-hiu„)\Ml 1 \ 



(8) 



For instance, the so-called broad band around the zero- 
phonon line at E — has a shape close to f(E)/E 3 at 
elevated temperatures (strictly speaking, the hole con- 
finement state participating in the transition has to be 
included as well). 



FIG. 2: Energy dependence of the phonon amplitudes for 
an InAs sphere of radius R embedded into GaAs material. 
The unrealistic situation without strain is compared with the 
strained one. 



III. RESULTS AND DISCUSSION 

In the following we present results for a spherical QD 
consisting of pure InAs embedded in infinitely extended 
GaAs. Under hydrostatic compression the coordinates 
change according to x'a = (1 — a)Xj. The deformation 
a can be determined from the static solution (lu = 0) 
of Eq. |21 Applying the boundary conditions from Eq. 0] 
accordingly leads to 



a = an 



C44,GaAs 



C44,InAs — C44,GaAs + jCnjnAs 



(9) 



with the lattice mismatch ag- For the present InAs/GaAs 
model system, it is a = (a ,i„A s - ao,GaAs)/ao,inA s = 
6.8% using the experimental values given in Tab. [I] The 
relaxed lattice constant a of the dot material (InAs), 
which is compressed hydrostatically, is therefore ai n A s = 
ao. InAs (1 — a) (last row of Tab.|U and indicated in Fig.Q] 
as dashed line.) 

Outside the dot, the material is compressed in radial 
direction but dilated laterally resulting in a zero net vol- 
ume change. Since we consider only hydrostatic com- 
pression the elastic properties are taken unchanged for 
the barrier material. 

We start by looking at the common approximation 
of using bulk phonons of the barrier material (index B) 
throughout. The solution of Eq. @ gives simply 



u {B) (r) = h{k B r) 



(10) 



with ji(x) being the spherical Bessel function of the first 
kind, and the wave nu mber i s related to the energy via 
E = hvisks (vi — \fC\xfp is the longitudinal sound 
velocity). Insertion into Eq. (|SJ) leads to 
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FIG. 3: Electron-acoustic phonon coupling function. Com- 
pared are (a) the unstrained combination of an InAs QD 
(R — 5nm) embedded in GaAs and the more realistic (b) 
strained situation. In each panel the assumption of bulk 
phonons, either GaAs or InAs, is compared to the full phonon 
calculation. The inset shows the electron band offset and the 
confinement charge density. 

In this approximation, the coupling function (divided by 
E 3 ) follows the (squared) Fourier transform of the elec- 
tron charge density, which decays on a scale of E « 
Uvib/R (R is the QD radius). 

For the full problem, the displacement inside the dot 
(index D) is getting a prefactor, 

u(r) = A(E)j 1 (k D r), (12) 

which has to be determined from solving Eq. El numer- 
ically. Although the integration in the matrix element 
extends into the barrier material as well, it is mainly 
the variation of A(E) which modifies the coupling func- 
tion. In order to implement the change in wave number 
(sound velocity) as well, we found it convenient to define 
a phonon amplitude as 

P(E) = A(E)—. (13) 

By definition, for bulk barrier phonons, P{E) = 1 which 
is the normalization used in Fig. [2J 

Whereas for a free standing sphere as well as for a 
rigidly clamped sphere, only discrete phonon energies are 
allowed, the QD embedded in an infinite elastic medium 
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FIG. 4: Electron-acoustic phonon coupling function: Com- 
pared are two different radii of the strained InAs QD embed- 
ded in GaAs. 
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[A] 


[A] 


GaAs 


0.067 


-7.17 





5.33 


5.65 


5.51 


InAs 


0.026 


-5.08 


-894 


5.66 


6.06 


6.05 


s-InAs 


0.026 


-5.08 


-208 


6.50 


5.79 





TABLE I: Parameters used in the calculation for the un- 
strained material 1 —. Given are also the theoretical lattice con- 
stants ao (DFT). Last row: Resulting values for the embedded 
strained InAs quantum dot (s-InAs). 



allows for a continuous manifold of phonon energies. 
However, the oscillating behavior of the phonon ampli- 
tude shown in Fig. |3 resembles the eigenmodes of an 
isolated sphere. A specific feature is the reduction of 
the phonon amplitude at small energies: Including the 
strain, its overall strength is reduced even further. Dou- 
bling the size of the embedded dot from R = 5 nm to 
R = 10 nm leads to halve the energetic spacing between 
maxima, while the overall amplitude remains the same. 
This can be understood by noting that within the ap- 
plied approximation, neither the static strain field nor 
the material density depends on the size of the QD. 

The calculation of the electron confinement wave func- 
tion is done within the effective mass approximation in- 
cluding mass discontinuities. All relevant input variables 
are taken from existing literatureii and given in Tab.[I] 
The compressive strain in the InAs dot reduces the con- 
duction band offset (CBO) significantly. The resulting 
electron confinement potential is plotted in the insets of 
Fig- El together with the charge density in the lowest con- 
finement state, r 2 ifii(r). In the strained case (Fig. El 3 ), 
a larger tunneling of the wave function into the bar- 
rier can be seen. However, this change in the electronic 
wave function affects the electron-phonon coupling only 
marginal. 

The coupling function Eq. (JSJ) of electrons with acoustic 
phonons is shown in Fig. El comparing the strain free case 
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(a) with the improved description including the strain 
present in the nanostructure (b). Within each panel 
the approximation of bulk phonons using either barrier 
(GaAs) or dot (InAs) material parameter are displayed, 
too. Compared to 'dot-bulk phonons', a strong reduc- 
tion is found, whereas 'barrier-bulk phonons' are much 
closer to the full solution, although missing the existing 
structure in the coupling function. 

The compressive strain results in an increased mass 
density (see Tab. Q), increased elastic constants in the 
dot and to a reduction of the coupling strength divu(r). 
The shape of the coupling function, however, is not much 
changed when including strain since the maximum posi- 
tions in the phonon amplitude are only slightly shifted. 
The coupling is reduced further by increasing the size of 
the QD, see Fig. Here, the coupling function is even 
more structured (sharp drop at around 1 meV) which 
has its origin in the oscillatory variation of the phonon 
amplitude discussed before. 

IV. CONCLUSIONS 

Including the details of the phonon characteristics in 
semiconductor nanostructures is essential for a quantita- 



tive description of the electron-acoustic phonon coupling. 
Material parameters, i.e. mass density and elastic con- 
stants, are changed due to intrinsic strain in pseudomor- 
phically grown heterostructures and are therefore spa- 
tially dependent. Within continuum elasticity this can 
be treated by introducing effective material parameters. 
There is a continuous energy spectrum for phonons of 
nanostructures embedded in an infinite elastic medium. 
Compared to the assumption of bulk phonons, an over- 
all reduction of the electron-acoustic phonon coupling is 
found due to the nanostructure as well as due to strain. 
Resonance-like features in the phonon modes lead to 
changes in the line shape of the energy dependent cou- 
pling function which are more pronounced for mid-size 
than for small dots. 
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